Self-stabilization of high frequency oscillations in semiconductor superlattices by 

time-delay autosynchronization 
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We present a novel scheme to stabilize high-frequency domain oscillations in semiconductor super- 
lattices by a time-delayed feedback loop. Applying concepts from chaos control theory we propose 
to control the spatio-temporal dynamics of fronts of accumulation and depletion layers which are 
generated at the emitter and may collide and annihilate during their transit, and thereby suppress 
chaos. The proposed method only requires the feedback of internal global electrical variables, viz 
current and voltage, which makes the practical implementation very easy. 
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Semiconductor superlattices have been demon- 
strated to give rise to self-sustained current oscillations 
ranging from several hundred MHz 0, 0, 3 to 150 GHz 
at room temperature [a. Various mechanisms with 
[E H H or without |10| the involvement of propa- 
gating field domains have been discusssed. In any case, 
a superlattice constitutes a highly nonlinear system > 
and instabilities are likely to occur. Indeed, chaotic sce- 
narios have been found experimentally 0, IT3. Ibi and 
described theoretically in periodically driven |15| as well 
as in undriven systems [lfj. For a reliable operation of 
a superlattice as an ultra-high frequency oscillator such 
unpredictable and irregular conditions should be avoided. 
In principle, synchronization of oscillations in a super- 
lattice by an external signal could be exploited to 
achieve a desired periodic behavior. However, in reality, 
the control of the forcing frequency in the ultrahigh range 
presents substantial technical problems. 

Here, we propose a simple self-stabilizing scheme that 
is especially suitable for semiconductor devices like super- 
lattices. It uses a profound concept of chaos control from 
nonlinear dynamics and chaos theory. Within this ap- 
proach, an intrinsically unstable time-periodic motion is 
stabilized using a simple feedback loop with a time delay 
18]. This type of control needs only small control forces 
initially, and they vanish once control has been achieved. 
A sound advantage is that the oscillation mode to be 
stabilized need not be known beforehand, in contrast to 
other chaos control schemes. Rather, a simple delay line 
leads to autosynchronization of the system. Methods of 
nonlinear control theory 0, [2(j have been usefully ap- 
plied to real world problems in various areas of physics, 
chemistry and biology [2l l I2I I2I I2I I2 I I2I I27I VM. |29 | . 
but no use has been made of this in the field of semicon- 
ductor self-oscillators. 
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Control methods can be either local or global 30] . Lo- 
cal methods require our ability to measure, and apply 
forcing directly to, the spatially resolved state variables 
of the system under study. However, in nanotechnol- 
ogy such variables, being e.g. electron densities in some 
quantum wells, are not easily accessible, and thus local 
methods cannot be applied. Unlike those, global methods 
require access only to some macroscopic variable(s) char- 
acterizing some integral output of the system. Such out- 
put can generally be reliably measured, and thus global 
methods seem to be the only option for the control of 
devices like superlattices. However, as we will show be- 
low, they are not straightforwardly applicable to nano- 
systems whose structure is spatially discrete. In this Let- 
ter we present a general approach to self-stabilization 
of irregular oscillations in semiconductor devices based 
upon essentially discrete quantum structures. 

We consider a model for nonlinear electronic transport 
in semiconductor superlattices that yields complex and 
chaotic dynamic behavior under fixed time-independent 
external voltage in a regime where self-sustained dipole 
waves jU are spontaneously generated at the emitter. 
Those dipole waves are associated with traveling field 
domains, and consist of electron accumulation and de- 
pletion fronts that in general travel at different velocities 
and may merge and annihilate. Such moving fronts are 
widespread in nonlinear, spatially extended systems, and 
similar chaotic front patterns occur in many other sys- 
tems, e.g., spatially continuous models describing bulk 
impurity impact ionization breakdown in semiconductors 
|33] or globally coupled heterogeneous catalytic reactions 
|33| . Thus the time-delay autosynchronization method 
proposed in this work could be readily applied to stabi- 
lize similar space-time patterns in a variety of systems. 

Our model of a superlattice is based on sequential tun- 
neling of electrons |31| . In the framework of this model 
the quantum wells are assumed to be only weakly cou- 
pled, and electrons are localized at these wells. The 
tunneling rate to the next well is lower than the typ- 
ical relaxation rate between the different energy levels 
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within one well. The electrons within one well are then 
in quasi-equilibrium and transport through the barrier 
is incoherent. The resulting tunneling current density 
Jm^m+i(F m , n m , n m+ i) from well m to well m+ 1 de- 
pends only on the electric field F m between both wells 
and the electron densities n m and n m +i in the respec- 
tive wells (in units of cm~ 2 ). A typical dependence of 
Jm-tm+i on the electric field between two consecutive 
wells is iV-shaped and exhibits a pronounced regime of 
negative differential conductivity. 

The rate of variation of electron density in well m is 
governed by the continuity equation 

e ~d~T = Jm ^ 1 ^ m ~ J™->rn+\ for to = 1, ... TV (1) 

and Gauss's law 

e r e (F m -F m -i)=e(n m -N D ) for to = 1, . . .N, (2) 

where N is the number of wells in the superlattice, e r 
and eo are the relative and absolute permittivities, e < 
is the electron charge, Nd is the donor density, and 
Fo and Fjf are the fields at the emitter and collector 
barrier, respectively. The total applied voltage U be- 
tween emitter and collector imposes a global constraint 
U = — X)m=o F m d, where d is the superlattice period. 
This, together with J2J), allows us to eliminate the field 
variables F m (ni, ...,tin,U) from the dynamic equations. 

At the contacts Ohmic boundary conditions Jo^i = 
ctFo, Jn->n+i — crFpfriN /Njj are chosen, where a is the 
Ohmic contact conductivity, and the factor n^/Njj is 
introduced in order to avoid negative electron densities 
at the collector. The value of a essentially determines 
the oscillation mode 

If the contact conductivity a is chosen appropriately, 
electron accumulation and depletion fronts are generated 
at the emitter. Those fronts form a traveling high field 
domain, with leading electron depletion front and trailing 
accumulation front. This leads to self-generated current 
oscillations. A fixed voltage U imposes a constraint on 
the lengths of the high-field domains and thus on the 
front velocities. If N a accumulations fronts and Nd de- 
pletion fronts are present, the respective front velocities 
v a and Vd must obey Vd/v a — N a /Nd- If the accumula- 
tion and depletion fronts have different velocities, they 
may collide in pairs and annihilate. At certain combi- 
nations of contact conductivity a and voltage U, chaotic 
motion arises, when the annihilation of fronts of opposite 
polarity occurs at irregular positions within the superlat- 
tice [l6|. The inset of Fig. ^ shows the plane of a and U, 
where regions with distinct regimes are marked by differ- 
ent shading. Black regions are those where chaotic be- 
havior has been found. As a computationally convenient 
criterion for chaos we have used the rapid decay of the 
autocorrelation function estimated from ri2o(t). Chaotic 
regimes are found at low contact conductivity and low 
voltages, where dipole oscillations with leading accumu- 
lation and trailing depletion fronts occur, and at higher 




FIG. 1: One-parameter bifurcation diagram: Time differences 
between consecutive maxima of the electron density in well 
no. 20 vs voltage at a = 0.5 Q~ 1 m^ 1 . Time series of 
length 600ns have been used for each value of the voltage. 
The inset shows a two-parameter bifurcation diagram: black 
squares denote chaotic oscillations, light shading indicate pe- 
riodic oscillations, and the white region shows the absence 
of oscillations. Simulation of an TV = 100 superlattice with 
Alo.3Gao.7As barriers of width b = 5nm and GaAs quantum 
wells of width w = 8nm, doping density Nd = 1.0 x 10 11 cm -2 
and scattering induced broadening F = 8meV at T = 20K. 



contact conductivity and higher voltage, where the role 
of accumulation and depletion fronts is interchanged. In 
Fig. ^ a one-parameter bifurcation diagram is given, ob- 
tained by plotting the time differences At between two 
consecutive maxima of the electron density in a specified 
well. The value of a is 0.5 f2 — 1 ni _1 , and U is changed. 
Chaotic bands and periodic windows can be clearly seen. 

The transition from periodic to chaotic oscillations is 
enlightened by considering the space-time plot for the 
evolution of the electron densities (Fig. At U = 

1.15V chaotic front patterns with irregular sequences of 
annihilation of front pairs at varying positions within the 
superlattice occur. We have calculated the largest Lya- 
pounov exponent as 1.1 x 10 9 s _1 which is a clear indica- 
tion of chaos. 

We shall now introduce a time delayed feedback loop 
to control the chaotic front motion and stabilize a pe- 
riodic oscillation mode which is inherent in the chaotic 
attractor. In the extended time-delay autosynchroniza- 
tion scheme first suggested by Socolar et al |34|, multi- 
ple time delays are used to improve the control perfor- 
mance. Analytical insight into those schemes has been 
gained only recently [35l l36l l3l| . and various ways of 
coupling of the control force, includin g lo cal and global 
schemes, have been compared [30L I38ll39ll4fj| . Whereas 
local coupling schemes usually lead to efficient control in 
a large control domain, they are not easily implemented 
in real systems since local, spatially resolved measure- 
ments are necessary. Therefore, here we propose a much 
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FIG. 2: Control of chaotic front dynamics by extended time- 
delay autosynchronization. a) Space-time plot of the uncon- 
trolled charge density, and current density J vs. time, b) 
Same with global voltage control with exponentially weighted 
current density (denoted by the black curve). Parameters as 
in Fig. □ U = 1.15 V, r = 2.29 ns, K = 3 x 10~ 6 , R = 0.2, 
a = 10 9 s~ 1 . Light and dark regions denote electron accu- 
mulation and depletion fronts in the space-time plots of the 
charge densities, respectively. 



simpler global scheme. In our problem, as a global out- 
put signal that is coupled back in the feedback loop, it 
is natural to use the total current density J defined as 

follows: J = TvTT ELo J"m->m+i H3- For tne uncon- 
trolled chaotic oscillations, J is given in Fig.|2K by grey, 
showing irregular spikes at those times when two fronts 
annihilate. Note that the grey current time trace is mod- 
ulated by fast small-amplitude oscillations (due to well- 
to-well hopping of depletion and accumulation fronts in 
our discrete model) which are not resolved in the plot. 
However, as the variable J is fed back to the system for 
the purposes of control, these high-frequency oscillations 
render the control loop unstable. They need to be filtered 
out by using e.g. the following low-pass filter: 



J{t) = a J(t')e"" (t ^'W, 



(3) 



with a cut-off frequency a. 

The multiple time delays of the resulting signal J 
(Fig. |21 black curve) are then used to modulate the volt- 
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FIG. 3: Control domain for global voltage control with ex- 
ponentially weighted current density. Full circles denote suc- 
cessful control, small dots denote no control. Parameters as 
in Fig. H 



age U across the superlattice: 

U = U + U c (t) (4) 
U c (t) = -K(l{t)-J(t-r)) +RU c (t-r) 

CO 

= -K R" (J(t - vr) - lit - + 1)t)) 



where C/o is a time-independent external bias, and U c is 
the control voltage. K is the amplitude of the control 
force, t is the delay time, and R is a memory parame- 
ter. Such a global control scheme is easy to implement 
experimentally. It is non-invasive in the sense that the 
control force vanishes when the target state of period r 
has been reached. This target state is an unstable pe- 
riodic orbit of the uncontrolled system. The period t 
can be determined by observing the resonance-like be- 
havior of the mean control force versus r. The result of 
the control is shown in Fig. \Bp. The front dynamics ex- 
hibits annihilation of front pairs at fixed positions within 
the superlattice, and stable periodic oscillations of the 
current are obtained. In Fig. [3] the control domain is 
depicted in the parameter plane of R and K. A typical 
horn-like control domain similar to the ones known from 
other coupling schemes [3(| is found. Typically, the left- 
hand control boundary corresponds to a period-doubling 
bifurcation, while the right-hand boundary is associated 
with a Hopf bifurcation. These findings show that our 
control scheme is robust. 

To conclude, we have demonstrated that time-delay 
autosynchronization represents a convenient and simple 
scheme for the self-stabilization of high-frequency current 
oscillations due to moving domains in superlattices. This 
approach lacks the drawback of synchronization by an ex- 
ternal ultrahigh-frequency forcing, since it requires noth- 
ing but delaying of the global electrical system output 
by the specified time lag. The proposed low-pass filter- 
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ing of the output signal presents a solution of the prob- 
lem one necessarily encounters when trying to control a 
nano-system with a crucially discrete quantum structure 
leading to superimposed fast well-to-well hopping oscil- 
lations in our case. 
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